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Abstract 

Discussions on a form of a frequency spectrum of wind-driven sea waves just above the 
spectral maximum continue during the last three decades. In 1958 Philhps made a conjecture 
that wave breaking is the main mechanism responsible for the spectrum formation. That leads 
to the spectrum decay u}~^, where u is the frequency of the waves. There is a contradiction 
between the numerous experimental data and this spectrum. The experiments show decay 
There are two general ways of the explanation of this phenomenon. The first one 
(proposed by Banner (1990)) takes into account the Doppler effect due to surface circular 
currents generated by long waves in the Phillips model. The second approach ascends to the 
work by Zakharov and Filonenko (1968). It is based on four-wave interactions in the kinetic 
equation and gives good agreement with the experimental data. 

In this article the contribution to the Phillips model due to the Doppler effect was examined. 
We considered the cases of isotropic and strongly anisotropic spectra. Both analytical results 
and numerical computations give us an evidence to state that the Doppler effect cannot explain 
the divergence of Phillips model with experimental data. 
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1 Introduction 



There is a long discussion in the scientific world about the form of wind generated sea wave spectrum in 
the "universal" region of frequencies uj > ujp, where ojp is the frequency of the spectral maximum. 

There are two main points of view. The first one proposed by Phillips in 1958 gives us the Philhps' 
spectrum in universal region 

Fi.) . ^, (1) 

where a is dimensionless constant, which according to experimental data appears to be rather small a ~ 
0.01. The second point of view introduced by Zakharov and Filonenko (1966) gives us in universal region 
following spectrum 

. (2) 

where /3 - is dimensionless constant, u - is wind velocity, e ~ ^ " is the ration between ocean water and 
atmosphere densities. 

The and (Q) formulas are based upon completely opposite propositions. The Phillips' spectrum (jl]) 
takes place if the spectrum at the high frequency region is mainly determined by appearance of sharp crests 
and wave breaking. The level of nonlinearity depending on steepness /i = fca in this case at least at the 
level of the wave peaks is extremely high - /i ~ 0.44. 

On the other side spectrum (||) is based on presumption of small mean steepness - /i <C 1. This proposi- 
tion is in great concordance with observed data. Usually /i < 0.1. 

Following this ideology the wind generated waves' ensemble is described by Hassebnan kinetic equation 
(Hasselman 1962,1963) 

^^Snl+P++p-. (3) 

Zakharov-Filonenko spectrum appears as exact solution of the particular case of Eq. ^ 

Snl = 0, (4) 

this is a classical case of Kolmogorov-Zakharov (KZ) spectra. This weakly turbulent theory is far advanced 
both analytically and numerically (Zakharov and Filonenko 1967, 1966; Craik 1987; Zakharov at al. 1992; 
Pushkarev and Zakharov 1996, 2000). 

A lot of experimental observations (Forristall, 1981; Kahma, 1981; Donelan, 1985, 1998, Hansen at al. 
1990) beginning from famous work of Toba (1973) gives us evidences that in most cases at moderately high 
frequencies ^ ^ 10 the Zakharov-Filonenko spectrum (^ is well detected. Sometimes in the region of 
extremely high frequencies ^ ^ 10 the Phillips spectrum ([l|) takes place. 

"Universality" of spectrum o;"^ became obvious in the middle of eighties. In 1986 O. Phillips, "father" of 
spectrum published an article with confession that spectrum (|l|) may takes place only at high frequency 
region and is not universal. 
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In the region of moderate frequencies the cj~* spectrum is dominating with confidence. From this fact 
the question of physical interpretation of spectrum F{uj) ~ uj^'^ becomes very important. The origin of this 
spectrum is rather obvious from our point of view. This spectrum is just weakly turbulent KZ-spectrum. 
There are another explanations, though. 

One of them was proposed by Phillips in 1986. He made supposition that spectrum u!~'^ is a result of 
balance of all three terms of Eq.(^ 



This affirmation unlikely can be treated as a serious theory, because the expression of is completely 
unknown. There are only very rough empirical formulas. However in that paper this expression was taken 
in the form giving us the balance in Eq. (||). But there is no any ground for this assumption. 

Another explanation was proposed by Banner (1990) and developed by Donelan (1998). They consider 
the uj~'^ spectrum as an artifact of the observation of time histories at a point brought about by Doppler 
shifting of the short waves riding on the orbital currents of the long waves. Let us suppose that spatial 
spectrum of the short waves $(fc) is known. The problem is to derive the frequency spectrum F{uj) using 
only ^{k). Normally we have to assume that waves obey linear dispersion equation. In this case the 
following simple expression takes place 



However in the presence of long waves the dispersion law of the short waves modifies due to Doppler effect 
at the orbital velocity's field v of the long wave component 



where v is the chaotic quasi-periodic function of time with the zero mean value. In this case the frequency 
spectrum F{uj) becomes wider. 

The present paper is devoted to the quantitative theory of this without any doubt very important phe- 
nomenon. We consider it in a very simplified way, supposing that v is sinusoidal function of time, i.e. 
assuming that the long wave is strictly monochromatic. This supposition is not too unreal. Undoubtedly 
this gives us the correct order of magnitude. All the more the result of this paper is quite clearly under- 
standable. 

In Sect. 2 we introduced the general notations (following Phillips). In Sect. 3 the general formula 
for spectrum (taking into account Doppler shifting) is derived. It is convenient to introduce a new set of 
dimensionless variables, this is considered in Sect. 4. In Sect. 5 we develop an important case of dispersion 
= gk ("deep water") where g is the gravity acceleration, k is the wave number equal to 27r/A. Section 6 
is devoted to numerical results. The conclusion is placed in Sect. 7. 

At the case of sufficiently small steepness /i ~ 0.05, the Doppler effect appears to be too weak to explain 
the clearly observed difference of ut^^ spectrum and experimental results, at least in the region of universal 
frequencies oj/ujp ~ 5 — 10. This fact makes too detailed inspection redundant. 



Snl ~p'^+p 



(5) 




L0{k) 



^W + (kv), 
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However, it is worth to say, that at higher frequencies the Doppler shifting can play very important role. 
This fact cannot influence on interpretation of experimental data at moderate frequencies though. At our 
standpoint the only conclusive interpretation of observed spectrum is treating it as weak turbulent 
Kolmogorov-Zakharov spectrum. It is worth to add that many direct measurements of spatial spectrum 
^{k) =< l^fcp > by photographing of the ocean surface and completely excluding the Doppler effect 
influence give us in the region of moderately high wave numbers the expectable result 

This fact is in great accordance with frequency spectrum F{ijj) ~ o;^** at the absence of Doppler shifting. 



2 Formulation of the problem and notations 

Following the notations of Phillips, the wave spectrum for homogeneous stationary wave field can be intro- 
duced as follows 

+00 

X{k, Lu) = (27r)-3 J J p(^r, t) exp[-i(kr - Lut)]drdt, (6) 



where p(r,t) = ^(x, io)C(^ + i"? ^0 + t) is the covariance of the surface displacement ^(x, t), r is the 



spatial separation vector, t is the time separation, k — (fci, fc2) — (fc, 6) is the wavenumber vector (the 
second representation is just the same in polar coordinates) and u; is the radian frequency. Obviously it can 
be rewritten in terms of Fourier transforms of corresponding functions 

(^(k,c^)r(k',c^')> = X(k,c^)(5k-k'<5^-^', (7) 

X{--k,-uj)^X{k,uj). 
It is convenient to introduce additional spectra which are the reduced forms of Eq. 

+CX) 

$(k) = 2 J X{k,uj)du;, (8) 



2-7T +00 

x{k,uj)kdkde. (9) 







In general case it is impossible to express F{uj) in terms of $(k). But in the case of small wave amplitude 
we can use the linearized equation 

Xik, w) = $(k) {Siuj -ujk) + 5{uj + ujk)) , (10) 
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where uJk — uj{k) is the dispersion law and S is the Dirac deha-function. Now one can rewrite Eq. in 
the following form 

F{uj) = 2j J ^{k,9){S{uj-ujk) + S{uj + ujk))kdkd0. (11) 


Let us remember that uit > 

+00 27r 

F{uj) = 2 I I <^{k,0)S{uj-ujk)d9kdk. (12) 





Using the rules of integration of (5-function one can obtain 

27r 



(13) 



k = k(Lu), 

where k{uj) is the inverse function for aj(fc). 
In the deep water case oj — \fgk 



g '''^''>~ 2y k~ 2u: 



27r 



^M = ^y (14) 



For an isotropic spectrum (<i>(/e, &) — $(fc)) 

Fi.) = . (15) 

For a strongly anisotropic spectrum ($(fc, 0) = ^(k)S(d — 6*0)) 

Fi.)^'44-\- (16) 



5^ V 5 

If the spectrum is determined by discontinuities of spatial derivative (wedges) caused by wave breaking 
according to Phillips (1958), we have 

^-^''fW-^- (17) 



In the isotropic case substitution Eq. ( 17) to Eq. (|15|) gives 



Fiuj) = Svra^. (18) 



Similar calculation for Eq. ( [T^ gives us the following result 



2 

FH=4a^. (19) 



Here a is the so-called Phillips' constant. 

If the wedges appear on each wave crest, a should not be too small (a ~ 0.1). Any attempts to fit Eq. 
dl^ or Eq. ( [19I ) to experimental facts give for a at least one order of magnitude less. 
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3 General formula for spectrum with the account for the Doppler effect 



Let us assume that we study waves on the background of periodic current v(t) beneath the surface. Taking 
into account the Doppler effect in the dispersion equation one can get the following substitution 



w(k) — >cj(k)+kv(t). 
Suppose that this current is described by a one-dimensional periodic function 

v{t) = V COSLUot. 

This allows us to rewrite our substitution in the following form 

LUk — > uj{k) + kv co8{ujot) cos ■0, 

where tp is the angle formed by k and v vectors. Without loss of generality one can take — 9. Using this 
expression one can write mean value of expression as 

T +00 27r 

F{uj) = J J J ^(.k, 9)6{uj - oj{k) - kv cos(woO cos e)dekdkdt, (20) 



where T is the period equal to 2tt/u!o. It is useful to take the average with respect to t as the first step. Let 
us denote this factor to the part of the integral as M 

T 

1 f 

M~ — f S(lj — kvcosoJotcosO ^ uJk)d(uJot). 

TT J 


The (5-function in this expression gives us the equation 

uj — ujk = kv cos 9 cos(a;oio), (21) 
This equation have the following roots 



LUotg — ± arccos 



kv cos 9 ' 



when the following inequality satisfied 



U) -U)k 



kv cos 9 



< 1. 



(22) 



(23) 



M 



TT\kv cos 9 sin {woto ) \ 
Using (|2|) and the last expression one can get the multiplier 

,x 2 2 

M = ■ 



kv cos 0a/ 1 



k'-^v'-^ cos-^ 



TTy/k^v"^ COS^ 9 — (W — Wfc)^ 



(24) 
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Using this expression one can get 

+ 00 27r 

It is necessary to understand that the limits of integration have to be chosen in order to satisfy inequality 



4 Convenient set of dimensionless variables 

It is convenient to introduce a new variable and a constant 

^^t}^^dk = ^dCX^-. (26) 
9 9 

One can see that C is a dimensionless analogue of k, and A is just the same for ui. In these variables the 
spectrum Eq. ( p^ takes the following relatively simple form 

fw-z^bT[^M5^, ill) 



^cos^e- (^^)' 



where B is the constant which appears as a result of the substitution k ^ in $(fc). For the Phillips 
spectrum case Eq. (|l7]) this constant takes the form B — av^ / g'^. 

The physical meaning of the constant A becomes more obvious if we transform it to the following form 

A^^^, (28) 

LOp Cp 

where usp — g jcp is the peak-frequency (this is the only frequency-like scale in the model) and Cp is the 
phase velocity. The ratio e = vjcp according to Banner (1990) is 0.05 in magnitude. 

For the case of strongly anisotropic spectrum ($(fc, ff) = ^{k)S{d)) one can easily obtain the following 
result by a simple substitution in Eq. (El) 



+00 

FiX)^'-^B f^M^. (29) 



Wl- 



If the spectrum is isotropic (^{k, 0) = <&(A;)) things become a little more complex. Following Eq. ( ^7| ) 
we have 

+ 00 7r/2 

F(A) = -4i? / $(C)dC / , (30) 
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An analysis of the 6'-part gives 

7r/2 



-l / cos^ f — 



, , 2 

2Q _ A-Afc 



c 



1 



K\,ll-[^)'y,K{^y (31) 

where K is the complete elliptic integral of the first kind. Using the inequality ( p3[ ) after substitution Eq. 
( pi] ) into Eq. ( pO[ ) one can obtain 

Expressions ( p9| ) and ( ^2| ) are valid for an arbitrary dispersion relation Aj, = uo{k)v/g. 
5 Deep water case 

For further study of the isotropic spectrum (this case is more complex so let us do it first and then perform 
the same calculations for the anisotropic case) we have to introduce the form of the dispersion law. We 
think that the case of the deep water surface is the most important. We have to use explicit expression of 
Wfc, in the case of a deep water to^ — y/gk or A/t — 



(33) 







Now we have to take into account condition ( P3[ ) 

C'- (a- VC)'>0. (34) 

Solving this inequality gives us 



yi,2 = - ± W - - A, 

2 V 4 (35) 
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where y — y/^, y > 0, so we can state that interesting domains are 

when A < — , 



F(A) 



dC + 



c 



dC 



(36) 



when A > 



F{X) 



8 9 



+00 



B 



(37) 



Let us explore special case, when $(C) is a quickly going down to zero function, like Phillips spectrum 
= In this case we can take into account only first part of ( at A ^ 1/4. Calculation of the 

spectrum decomposition is situated in Appendix ^ and ( A7) gives us the following result 

1 



F(A) = Fo(l + ^A2 



+ ...), 



where 



v"" 1 
7^' 



(38) 



(39) 



this is just the Phillips spectrum Eq. (y_8|) in absence of Doppler effect for isotropic distribution. 

For the one-dkectional spectra we can use exactly the same approach, but all calculations become simpler. 



The starting domains of integrations (which correspond to Eq. 
1 



and Eq. (b7b) are the following 



when A < 



F(A) 



2 9 



B 



dC 



c 



+ 00 



HC)dC 



-A 



1 - 



A-n/C 

C 



(40) 
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when A > -, 
4 

-\-oo 

FiX)^i<LB I (41, 

The complete analysis of this case is situated in Appendix ^ and (A9) gives us the following formula 

F{X) ^ Foil + ^X^ + ...), (42) 

where 

1 

Fo = 4a— — , (43) 

this is just the Phillips spectrum Eq. (|9]) in absence of Doppler effect for strongly anisotropic distribution. 

Let us try to evaluate the range of applicability of the obtained results. As it was postulated in Eq. ( p8| ) 
we can state 

LO V 1 

A = < -. 

luq Cp 4 

Following Banner (1990) we can take e = v/cp = 0.05. Using this one can obtain the following 
inequality 

u) < 5wo- (44) 

It means that we cannot consider the uj~^ spectrum as an artifact of the observation of time histories at a 
point brought about by Doppler shifting of the short waves riding on the orbital currents of the long waves. 



6 Numerical results 

The integrals Eq. Eq. Eq. (|^ and Eq. can be easily calculated numerically. It is not very 
illustrative to view the explicit form of the spectrum, so let us represent it as the ratio i^(A)/Fo(A). The 
expressions Eq. (^6|), Eq. ( p^ and Eq. (|7|), Eq. (^ij) give us the curves shown in Fig. |l[ 

It is obvious that the form of the first part of the decomposition is not linear and it means that we cannot 
transform the u!~^ to u''^. The most illustrative form to estimate the power of the ui in F{X) is the logarithm 
scale. The results are represented in Fig. § 

As one can easily see, there is almost no difference betwin uj~^ line and calculated curve. Thus in this 
range (A < 1 /4) we can state that the does not give Doppler a desirable result. We can see that there are no 
prerequisites for transformation Phillips spectrum (u;^^) to experimentally observed tj^^. 



7 Conclusion 

The presence of at least two different approaches on the wind-driven sea waves frequency spectra formation 
motivates us to examine the Phillips model in the presence of the Doppler effect. The pure Phillips model 
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gives the results essentially different with the experimental data. Following Banner conjecture we obtained 
the general formula for an arbitrary frequency spectrum in the presence of the Doppler effect. Using this 
formula we have performed an analytical analysis of the influence of the Doppler frequency shifting on sur- 
face waves' spectra at the deep water and strongly anisotropic cases. The forms of spectra were calculated 
numerically for a wide range of frequencies. 

Analytical and numerical results give us an irrefutable proof of inconsistency of Doppler shifting mech- 
anism in fitting Phillips model with the experimental results. An analytical estimation shows that in the 
presence of Doppler effect its contribution is negligible. The numerical computation gives spectrum decay 
~ u!^^ at least up to frequency ~ 20uip. 

According to the present experimental results one has to recognize that PhUUps model of energy dissipa- 
tion does not give the result similar to experimental data even if the Doppler effect takes place. 



Appendix A Calculations of the spectra decompositions 

At first, a Uttle bit more complex case of isotropic spectrum. Let us perform the following substitution 

c ■ 

In this case in our domain of integration (— 1 < a; < -|-1) we can obtain the invert representation 



1 - 2\/l + 4a;A -|- 2xX 



(Al) 



(A2) 



and further 

(^(x) « ^ [2X^x'' - AX'x'' + WX^x^] =X^- 2X^x + hX^x\ 
For substitution dC, = C,'{x)dx we have to calculate the first derivative 



(A3) 



, 1 A VI + 4Aa; 



a;3 a;Vl + 4Ax 



_ A ~ _ J_ _ A n _ 2Aa; + QX^x^ - 20X'x^ + 

+ 70A*a;^l + +^ \l + 2Xx - 2X^x'^ + AX^x^- 

- lOX^x" + 28A5a;5] - A = 
■' X'^ 

= -2A^ (1 - 5Aa; + 2lX^x^) . (A4) 

$(C(a;)) = $(A2(1 - 2Aa; + SA^a;^)) « 

«$(A2)-5lp2A3a; + 

+ 51^5AV + qp4AV. (A5) 
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Now we can calculate the integral ( pq ) with Phillips spectrum (<I>(C) = 

+1 



F(X) = --^B2X^ / <^>(C(x))C (x)K(^/l - x^)dx, (A6) 

-1 

Taking into account several well known results 

+1 



2 



-1 
+1 



xK {\/l-x^]dx = 



+1 +1 



-1 

where is the complete elliptic integral of the second kind, one can obtain the final result 

f^(A) = 8^4^(1 + iA2 + ...). (A7) 



Calculations of the spectrum decomposition for strongly anisotropic case up to expression Eq. (46) are 
just the same. Integral ( |4C| ) takes the form 

F(A) = H^2A3 fiipm.., (A8) 

-1 

The final result is 

F(A) = 2a^^(l + iA2 + ...). (A9) 
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Figure Captions 



Fig. 1. Ratio of calculated spectra F to Phillips spectrum as a function of dimensionless variable A is shown. It is 
clearly that F starting to deviate from Phillips spectrum only at A ~ 0.4 

Fig. 2. Logarithm of calculated spectra F and Phillips spectrum as a function of logarithm of dimensionless variable 
A is represented. It is easy to see that calculated spectra are lines are exactly parallel to line corresponding to Phillips 
spectra. This allow us to conclude that the power of frequency is the same and is equal to -5 
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Figures 




Fig. 1. Ratio of calculated spectra F to Phillips spectram as a function of dimensionless variable A is shown. It is 
clearly that F starting to deviate from Phillips spectrum only at A ~ 0.4 
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Fig. 2. Logarithm of calculated spectra F and Phillips spectrum as a function of logarithm of dimensionless variable 
A is represented. It is easy to see that calculated spectra are lines are exactly parallel to line corresponding to Phillips 
spectra. This allow us to conclude that the power of frequency is the same and is equal to -5 
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